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and Thomas A. Byrdsong 


SUMMARY 


This paper presents the results of a preliminary investigation of 
the effect of nonlineair structural terms on the flutter of a two-degree - 
of -freedom system. The three types of nonlinearities Investigated were 
a flat spot, hysteresis, and a cubic spring. Calculations were made on 
an analog conputer. For one case, the flat spot, an experimental inves- 
tigation was also made and good correlation with theoiy was found. 

In general, it was found that the linear flutter speed did not 
change for small disturbance angles} however, for larger disturbance or 
input angles, the flutter speed itsually decreased. One exception was 
the cubic hard spring, for which a limited-amplitude flutter was fovmd 
to exist we3JL above the linear flutter speed. 


INTRODUCTION 


Ifetil fairly recently the problem of aircraft flutter has been 
treated analytically as a pinrely linear phencanenon. Potentially, how- 
ever, many sources of nonlinearities exist and their possible effects 
are receiving more and more attention. It is thought, for example, 

' that nany of the limited-anplitude oscillations are associated with 
nonlinearities . 

In the flutter problem three types of forces - the inertia, the 
danping, and the elastic - are Involved. Each of these forces may arise 
from either the aerodynamics of the problem or from the structure. Non- 
linearities associated with the aerodynamics may arise, for example, 
from the effects of thickness or flow separation or frau the variation in 
hinge mcanents with anplitude. Some aerodynamic effects of thickness 
have been studied analytically by VanTDyke (ref. l) . He has treated the 
case of two-dimensional, supersonic flow and has found, except near a 
Mach number of 1, that thickness effects for this case are moderate. 

GThe structural side of the problem may also give rise to nonlin- 
earltles in the inertial, damping, and elastic forces. Thp present 
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invest igat ion is concerned with, nonlinearities in the elastic forces^ and 
three types of nonlinear springs have been considered. Their effects on 
flutter have been examined in the investigation of a two-dimensional, two- 
degree -of -freedom system. A part of the investigation has been made on an 
electronic analog computer. The analog representation of the problem is 
discussed in the appendix. 
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SYMBOLS 

natTjral frequency in translation, radians/sec 

natural frequency in pitch for basic linear system, radians/sec 

mass moment of inertia about elastic axis per unit leng;th of 
span, slug-ft^/ft 

spring constant of linear system, ft-lb/radian 

moment, ft- lb 

mass per unit length, slugs/ft 

maas static -unbalance per unit length, referred to elastic 
axis, slvig-ft/ft 

flutter speed 

flutter speed of 3JLnear system ' ' 

coordinate of cen-ter of gravi-ty measured from leading edge, 
percent chord 

coordinate of elastic axis measured from leadi n g edge, per- 
cent chord 

angular free play, deg 

angle of rotation, deg unless Indlca-ted otherwise 

flvild densl-ty, sl\igs/cu ft 

nondlmensional distance of pitch axis from midchord measured 
in half-chords, positive for positions of pitch axis behind 
midchord, a = Sx^g, - 1 

wing semichord, ft 

veidiical displacement from equilibrium, ft 

time, sec 


Dot over qxiantity Indica-tes differentiation -with respect to time. 
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ITXPES OF NOHLINEAR SPKENGS IMVESTIGATED 


Nonlinearities in tbs elastic forces can arise in many possible 
ways. Three types of nonlinear springs have been considered, and 
their characteristics are shown as plots of the force or moment required 
to produce a given displacement. (See fig. 1.) 

Shown at the left of the figure is a flat spot. This may be con- 
sidered to represent the condition of free play in the hinge or linkage 
of a control system. 

At tbe right of the figure is one type of hysteresis. As force or 
moment is increased, displacement varies in a linear manner mtil a 
point is reached at which a junqp occurs, after which the system is again 
linear. On the return path, a corresponding jump occurs at a negative 
value of the force. A nonlinearity of this type may occur in the case 
of a control surface with free play if friction exists at some point in 
the linkage, or it may occur through the action of power servos. In 
the case of a wing, this type cJf nonlinearity might represent the effect 
of rivet slip. 

The third type of nonlinearity treated has been termed a ciibic 
spring, where the force exerted depends on the usual linear power of 
the displacement and, in auidition, on a term containing the cube of the 
displacement. This can be considered either as a hard spring, which 
becomes stiff er as displacement Increases, as shown by the solid line, 
or as a soft spring, which beccmies weaker as displacement increases, as 
shown by the dashed line. In the case of a control surface, a hard- 
spring effect might be associated with power controls. In the case of 
a structure, a hardening effect is found when a thin wing, or perhaps 
a propeller, is s\ibjected to increasing amplitudes of torsion. A soft- 
spring effect may be associated with panel buckling. 


EFFECTS OF NONLINEAR SPRINGS ON AIRCRAFT FLUTTER 


These three types of nonlinearities have been considered in the 
fltctter investigation of a two-dimensional, two -degree -of -freedom system, 
free to oscillate in translation and in pitch. The flutter character- 
istics of the system have been determined with an electronic analog 
computer, vising Inconrpressible, indicial air forces. GJhe. analog repre- 
sentation of the fundamental linear flutter equations was based on the 
work of Baird, Pines, and Winson (ref. 2 ) . For the flat spot, in 
addition to the analog results, experimental results have been obtained 
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in tlie Langley 2- "by It-foot fluirber research -tunnel. The physical 
characteristics of the systems -treated for each of the three -types of 
nonlinear springs are given in table I. 

In the flu-tter results to he discussed, one of the -variables to he 
enployed is an input angle 0. The significance of this Inpirt angle is 
given in the following discussion. 

In ob-taining results with the analog ccmputer for a partlctilar case, 
a given "value of velocl-ty is selected, the system is given an Initial 
displacement (for exanple, 0) , and i-bs response to this displacement 
is examined for stabill-ty. In geneiul, the results to follow will shew 
that the response of the system becomes unstable only when a sufficiently 
large initial displacement is inposed. In the plysical system, this 
initial displacement corresponds to the effect of a gust or to an abrupt 
movement of the control stick. 


Flutter Results With the Cubic Spring 

In figure 2, analog results are shown for the cubic spring in the 
torsional degree of freedom. Results are shown in the form of the input 
angle 0 required to induce flutter as a function of veloci-ty. For 
this case, the relation he-tween moment and displacement was arbitrarily 

assumed to be M = lll0.59 -t ^30^, where M is the moment in foot- 
pounds and 0 the torsional displacement in radians. The coefficient 
of the linear term in 0 represents the spring constant in foot-pounds 
per radian of the linear system. 

First, the flutter speed of the linear system -was investigated and 
was fo\md -to be around 27 O feet per second. Then, both hard a-nfl. soft 
cubic springs were considered. Flutter boundaries are shown by the solid 
curve for the soft spring and by the dashed Hup for the hard spring. 

In both cases, the flu-tter region lies to the rlg^t of the boundary. 

For the hard spring, the flutter boundary is a s-trai^t l-tuf* at the 
flutter speed of the linear sys-tem. The soft spring in this case had a 
des-tabilizlng effect in that flutter could be induced below the linear 
flutter speed by making the initial, displacement sufficiently large. 

For the system treated herein, however, the de-vlation from the linear 
flutter speed occurs only at fairly hi^ initial displacements. 

It is also of interest to consider the flu-tter au 5 )li-tudes associated 
with these results. For both the linear system and the soft spring, 
the fltrtter oscillation -was highly divergent at any veloci-ty above the 
fltrt-ter boimdary. With the hard spring in the sys-tem, however, the 
fl-utter amplj-tude is self limited. This limit arapll-tude is a function 
of veloci-ty, however, and increases as velocity is increased beyond the 
fl-utter boimdary. 
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In studying other configurations, cases have been found in which 
a hard spring can he destabilizing, and it appears that the effects 
produced by a cubic spring depend on the stiffnesses of the original 
linear system. Generally, flutter speed decreases as the bending- 
torsion frequency ratio approaches unity. If a ciJbic spring, whether 
hard or soft, tends to mke this ratio approach unity, it will probably 
be destabilizing. 


Flutter Results With the Flat ^ot 

The second type of nonlinearity to be considered is the flat spot 
and, as mentioned previously, both tJind-tunnel and analog studies were 
made. The main features of the model used in the wind-tunnel tests are 
shown in figure 3* The sketch represents a two-dimensional model which 
completely spans the test section and which is free to oscillate in 
pitch and translation. The pitching degree of freedom is provided by 
a bearlng- 8 t 5 )ported shaft located slightly behind the quarter chord. 

This system in turn is suspended between a pair of leaf springs on 
either side of the test section so that the entire mechanism is free 
to translate. 

In translation the system is linear. The nonlinearity was intro- 
duced in the torsional degree of freedom .in the manner shown in the 
detail sketch. A leaf spring is clamped to the end of the torsion axis 
and its free end extends tpward between two set screws. The gap between 
the screws can be closed completely to give a linear torsion spring or 
opened to provide any desired amount of angular free play, giving the 
spring characteristic shown to the right in figure 3» 

In figure 4, results obtained in the wind tunnel with this conflg- 
xiratlon are conpared with analog results. The ordinate of the figure 
is again the input angle 9 in degrees. The abscissa is a flutter- 
speed ratio where is "tbe flutter speed of the linear 

system, and it is immediately seen that the flutter speeds have been 
decreased. In the investigation the angular free play brp was kept 
constant at 0.5°* The solid curves and the regions labeled stable, mild 
flutter, and flutter refer to the analog results. The circles show wind- 
tunnel experimental points. 

In both the experiment and the analog results, a preload was incor- 
porated into the system. This prelosid was comparable in effect to a 
deflected tab and produced a moment which varied with the velocity. At 
equilibrium, therefore, the model was not centered in the flat spot but 
rather at some point on the linear arm of the diagram. Ihis means that, 
for small, displacements from equilibrium, the system moved only on this 
linear arm so that, at velocities below the linear flutter speed, the 
system was stable. 
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With displacements slightly greater than 0.5°^ ‘fche effect of the 
flat spot began to enter the picture and. a region of mild., limited- 
amplitude flutter vras encountered. Here the system oscillated about 
its equilibrium position across only one knee of the flat spot and did 
not traverse the full flat-spot width. 

As an example, at a value of the flutter -speed ratio of 0.6, the 
flutter amplitude in this region of mild flutter was about 1°. Such 
limited oscillations, when translated into the behavior of control 
surfaces in flight, might not be dangerous in themselves but could be 
significant from a fatigue standpoint. 

As initial inputs were increased still further, a much luore violent 
flxitter was encountered in which the model osc il lated conpletely throu^ 
the flat spot with large amplitude. 5Che experimental wind-tunnel points 
apply to this violent type of flutter and agreement with the analog 
results is fairly good. Although not shown in the figure, it should be 
noted that in the esperimenteil tests^ also, a region of mild flutter 
was encountered at Initial displacements conparable to those indicated 
by the analog. 

In obtaining these res\ilts, only a single value of the preload has 
been taken into account, and further studies in which the preload is 
varied would be desirable. The destabilizing effect of the flat spot 
in this case is probably associated with an effective reduction in tor- 
sional stiffness. In the absence of any preload, the torsional stiff- 
ness at equilibrim would be reduced to zero, and presumably then one 
would at least expect to encounter tmstable osclUations in the region 
of the flat spot. 


Fltctter Results With ^steresls 

If, in addition to the free play, some source of static friction 
exists, the hysteresis phenomenon occurs. As mentioned previously, a 
hysteresis of the type described may arise in the case of friction in 
the hinge or linkage of a control surface or, in the case of a wing, 
might be associated with rivet slip. 

Before considering the flutter results with hysteresis, seme purely 
structural characteristics of the system with hysteresis should be exam- 
ined. Shown in figirre 5 is the response of a system with hysteresis at 
zero airspeed. The upper and lower traces represent, respectively, the 
displacement of the system and the variation in moment after release 
from the initial displacement fraa equilibriimi. The indentations on 
the lower trace occur when the system passes through the flat spots of 
the hysteresis box. 
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The presence of this hysteresis in the system has two primary 
effects on the structural characteristics. One is the introduction of 
variable structural damping, shown by the lines drawn tangent to the 
upper trace. At high amplitudes, the rate of decay is relatively small. 
As the an 5 )litude of the oscillation approaches the height of the hys- 
teresis box, the dangling is considerably greater. i-/hen the oscillation 
is contained within the box, the structural damping in this illustration 
is zero. 

The second effect of the hysteresis on the structural properties 
of the system is not apparent in this figure but is the introduction of 
an effectively weaker spring. At high angilitudes the frequency of the 
oscillation is less than the frequency at low ainplitudes where the 
system is linear. 

Analog results for bending-torsion flutter with hysteresis in the 
torsional degree of freedom are shown in figure 6. Results are again 
presented in the form of the input angle 0 plotted against flutter 
velocity. Two valiies of 6^, the amount of free play in torsion or the 

width of the hysteresis box, were considered. These were 8^ of 0.2°, 
shown by the solid line, and 8 qi = 0.4°, shown by the dashed line. The 

height of the hysteresis box was held constant. In these results, the 
flutter region lies above or to the right of the lines. 

Consider first the results for 8^ of 0,2°, shown by the solid 

c\irve. For very small inputs, the system oscillates on a line through 
the center of the box with a linear spring constant, as shown by the 
dashed line in the spring diagram. For this reason, the flutter bound- 
ary at the highest velocity occurs at the flutter speed of the linear 
system. 

This boundary continues straight upward until an initial displace- 
ment of 0.6° is ingjosed. Actually an initial displacement greater than 
0.2° represents a displacement beyond the linear path within the hys- 
teresis box. At low amp litudes, since the presence of the hysteresis 
introduces a large amount of structural d a m ping, for inputs up to 0.6° 
the oscillation decays and the system oscillates at constant aog)lltude 
within the hysteresis box. For larger displacements at the linear 
flutter speed a.ud at velocities above the linear flutter speed, the 
oscillation diverged rapidly. 

As in the case of the flat -spot results, flutter could be induced 
at velocities below the linear flutter speed by niald.ng the initial 
displacements sufficiently large. At the larger displacements, the 
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hysteresis dandling is smaller and the presence of the free play in the 
system is producing an effectively weaker spring. Finally a velocity 
is reached helow which the -air forces will not sustain an oscillation 
regardless of the magnitude of the initial displacement. For all veloc- 
ities helow the linear flTxbter speed, the flutter is of limited amplitude. 

Results for the wider hysteresis hox exhibit the same trends. At 
velocities below the linear flutter speed, however, greater initial dis- 
placements are required to Initiate fltrtber, since the wider hysteresis 
box introduces more danping. At the same time, the increased amount of 
free play reduces the effective spring constant even more, and, the 
region of limited-anplitude flutter extends to a lower velocity. 


COHCLTJDING EMARKS 


Little is known of the exact nature and magnitude of the nonlinear- 
ities which exist in actual aircraft. The purpose herein has been to 
make a prelimlnaiy sttidy of the effects of some nonlinearities which might 
occur. The results shown are, of course, a function of the particular 
configurations treated and are perhaps most closely related to the prob- 
lem of control-surface flutter. They indicate, however, that in some 
cases nonlinear effects can be large and can be destabilizing . The 
resvilts indicate further that the stability of a nonlinear system can 
become a function of the magnitude of an external excitation; it is 
also indicated that, when a nonlinear system becomes ixostable, its 
flutter may become less violent and self -limited. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May S, 1955 • 
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APPENDIX A 

ANALOG REPRESENTATION OF NONLINEAR FIUTTER EQCIATIONS 


In this appendix the wiring diagrams for the various types' of non- 
linearities which have been treated are considered. In figure 7^ a 
diagram is given -vdiich shows the wiring for the analog representation of 
the equations of motion for the conqplete system with the dashed hox indi- 
cating the point at -vdiich the nonlinearities are introduced. Figures 6 , 

9, and 10 show, respectively, the wiring for the cubic spring, the flat 
spot, and the tiysteresis nonlinearity. The function 0^ shown in fig- 
ures 9 and 10 denotes the output of the nonlinear circuits. 

The form of the basic equations of motion used in the wiring is 
essentially that given by reference 2. In reference 2 the equations of 
motion are general in that they eire written in terms of coupled modes of 
vibration and provide for spanwise variations in structural properties 
and in deflections. In the present analysis a two-dimensional system 
oscillating in rigid modes of translation h and pitch 0 is considered. 
For this case the equations of motion can be written in the following 
foim; 


For the translational degree of freedom. 


-h = 




1 . 0 2 .. 1 , 0 
■^ 11 ^ ^ 12 ^ 




A^0(T)h(O) + A^0(T)h(O) + A^0 (t)0(O) + A^0(t)0(O)J (A1) 
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and for the pitch degree of freedom. 




22^ 


4i^ 


4^ 


4- 


Ai,« 

22 


+ A 


0 

22 


0 + 


5 5 

A21 / ‘ i '^1 + ^22 

J 0 

4i fa 


fa «’("-"i)®("-"i) ^"1 
fa ("-" 1 ) ^"1 


5 ^ j • 

A^0(T)h(O) + A23^{25(T)h(0) + 


A220(t) 9(O) + App0(T)e(O) 


22^ 


) 


+ 


+ 


(A2) 


where t = Vt/b is a nondimens ional time factor and is the varia- 
ble of integration. For the cases considered, the quantities appearing 
in these eqxiations are defined as 


2 2 
A^ = m + rtph 

A^ = 2rtph^ 


0 /tLvbxS 
^11 ” \ w / 


m 




1 X 

A^ = 2itph^(l - a) 


0 3 

= 2itph 

4 = ° 

= -atpb^ 

4 2 

A^ = -2rtph^ 


= -2.pb5(| - a) 
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The term 9 J(t) is the "lift-deficiency function" obtained from the 
approximate expression by Robert T. Jones (ref. 5) for Wagner s indicia! 
response fionction (ref. 4) and appears as 

0(t) = + a2S^^^ (A^) 


where 

a^^ = 0.165 = 0.535 

p3_ = -O.Olt-55 p2 = 

The nonlinearities were introduced in the pitch degree of freedom 
and were represented by circuits which correspond to the y>ring- 
characteristic diagrams in figure 1. The nonlinearities are associated 

with the first term of ^ 2 ' ^ case of the cubic spring (fig. 8) 

the torsional frequency aw appearing in A® was expressed in terms 

5 ^ 

of the moment function M = Kq 0 + C0 which was assumed for the spring- 
characteristic diagram of figure 2. 

Potentiometer settings for the diagram of the equations of motion 
in figure 7 are given in table II. 
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TABLE I.- CHARACTERISTICS OF SYSTEMS 


Parameter 

Systems 

Cubic 

1 

Flat spot 

^steresis 

Semlchord, ft 

percent chord 

Xgg, percent chord 

alxigs/ft 

slug-ft/ft 

I^^, sliig-ft^/ft 
03 ^, radians/sec 
radians/sec 
Q, sl\igs/cu ft 

0.5 

23.0 

25.6 

O.lJli -89 

0.003871 

0 . 00893 i«- 

9 ^t .25 

125.66 

0 .002378 

0.5 

29.7 

51.8 
0.2985 
0.00521 
0.01196 

58 . 9 '^ 

81.24 

0.00245 

0.5 

23.0 

25.6 

0.1489 

0.005871 

0.008934 

66.60 

125.66 

0.002578 
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TYPES OF NONLINEAR SPRINGS INVESTIGATED 
BDRCE 





CUBIC 


Figure 1 


FLUTTER WITH CUBIC SPRING IN TORSION 



Figure 2 
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NONLINEAR FLUTTER MODEL 



FLUTTER WITH FLAT SPOT IN TORSION 



Figure 4 



c 


MCA TN 3559 


1 ? 


DAMPING ASSOCIATED WITH HYSTERESIS 




Figure 5 


FLUTTER WITH HYSTERESIS IN TORSION 



Figure 6 










